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Abstract 

The diffusion equation which describes the evolution of the average one 
particle distribution function for an ensemble of strongly turbulent plasmas is 
derived. The diffusion tensor is a time integral of the autocorrelation tensor of 
the fluctuations as observed by particles moving along statistically distributed 
orbits. These orbits contain the effects of fluctuations and thus differ from those 
encountered in weak turbulence theory. Two statistical orbit effects quadratic 
in the strength of the fluctuations affect the magnitude of the diffusion: (a) modi- 
fication of the ensemble average orbits by the fluctuations and (b) statistical dis- 
persion in particle orbits about the average. The plasma trajectory equations 
are used to relate each to the diffusion tensor itself when the turbulence is elec- 
trostatic. The diffusion tensor is explicitly evaluated for a strongly turbulent 
unmagnetized plasma. 
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RESONANT DIFFUSION IN THE PRESENCE OF 


STRONG PLASMA TURBULENCE 


I. INTRODUCTION 

Plasma turbulence is described theoretically by statistical formalisms: 
one derives equations which advance in time the average properties of an en- 
semble of plasmas. The ensemble may be thought of as consisting of realiza- 
tions which differ only in the phases of the Fourier decomposed microfields at 
some arbitrary instant. Rigorous specification of the ensemble is, however, 
generally unessential to the formalisms. 

For each realization, the one particle distribution function f for each species 
satisfies the Vlasov Equation 


>£-I-o (1 > 

with F the force per unit mass on the element of plasma located at the phase 
space point x , v at time t. Denoting the average of a quantity over the ensemble 
by < > and its deviation from the average by 8 we obtain from Eq. (1) 

ft + t<f? = - <££■ !f^> (2) 

By subtracting Eq. (2) from Eq. (1) the equation for 8f, the fluctuating com- 
ponent of the distribution function in any single realization, is obtained 


1 



w 

At 


v-zSf f (<E>tSE)’- ft 




In the usual theory of weak turbulence 1 ' 9 one first discards terms quadratic 
in the fluctuations from Eq. (3) and then integrates the resulting linear equation 
for S f . Plugging this expression for 8 f into Eq. (2) and making standard assump- 
tions result in a diffusion equation for < f > 


2<f7 


t- + XE?' 


3 . tj . ?S£ 

V Jv ~ 2V 


(4) 


The diffusion tensor in the weak turbulence approximation is 


°° X 

P -- jdT- F(><, V, t) J’ff? y t-r]/ 


o 

where ** and v* are the characteristic trajectories of the linearized Eq. 3, 
solutions to the ordinary differential equations 


(5) 






( 6 ) 


with the boundary conditions x*(t) =x , v* (t) = v. Note that the orbit x* (t), 
v* (t) is a non- fluctuating one. 

By taking moments of Eq. (4) one obtains equations which describe the evolu- 
tion of the ensemble averaged macroscopic properties of the plasma. 

From Eq. (5) we see that the magnitude of the diffusion tensor element D 
depends on both the strength of the fluctuations and the effective width (as a 
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function of r) of the correlation tensor element 


Cl j -- (*,!', t) 

where SF. is the i th vector component. The net area under the C.. vs. t curve 
divided by I SF[ ms | 1 SF T ms | i s a measure of this correlation width. Except for 
particles with approximately the same velocity as the typical phase velocity of 
the waves which comprise the turbulence, the net area under the C vs. r curve 
is small (Either C is a rapidly damped or highly oscillating function of t.) be- 
cause the turbulence is convecting past such non- resonant particles. 

The correlation width and D. . are largest for resonant particles. 10 For 
resonant particles C has the effective width r c = (Aw*)" 1 , where Aw* is the 
frequency spread in the Fourier components of the turbulence as observed by 
the resonant particles. Thus r c is the correlation time of the fluctuations as 
apparent to the resonant particles moving as if the fluctuations were absent. 

(We shall define r c as such even when the fluctuating fields are large.) 

In a strongly turbulent plasma the kinetic equation for <f> is a diffusion 
equation. It is again resonant particles which are most strongly diffused. Our 
treatment in this paper is concerned exclusively with them. The diffusion tensor 
becomes significantly different from Eq. (5) once the fluctuations reach such a 
magnitude that in the time interval r they appreciably alter particle orbits. 

In our theory of strong turbulence the effect of fluctuations on particle orbits 
is twofold. Consider a particle which in each realization of the ensemble is at 
the same phase space point at time t. As time progresses <x*> , <v*> — the 
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phase space orbit of this particle averaged over the ensemble — develops a 
contribution quadratic in the fluctuations analogous to the frictional drag in a 
Fokker- Planck description of interacting particles. [This contribution occurs 
because in integrating x* = F, it is F (x*, v*, r) which enters. Although 
<$F (x,v,t)>= 0, <SF(x*, v*, r) > ^ 0.] There will also be a statistical 
spread in orbits about <x *>, the dispersion likewise being an effect quadratic 
in SF. 

Suppose the motion of the particle at time t is resonant. With advancing 
time <v*> changes and the particle gradually becomes non-resonant. Were the 
particle to become non-resonant in a time interval r l shorter than r c (We esti- 
mate Tj to be the time needed for the S F part of <x * > to grow to order 27r/Ak, 
the characteristic scale length of the turbulence.), then the width of the auto- 
correlation as viewed by the particle would be Tj rather than r c . (The wave 
number spectrum of the turbulence has width Ak and has typical wave number k Q .) 
The fluctuations may similarly act to make an initially non- resonant particle 
resonant. 

The spread in orbits also acts to modify (in the ensemble average) the auto- 
correlation time of the fluctuations as viewed by the particle. P [x*(r), v* (r), 
r; x, v, t] , the probability of finding the particle in any single realization at 
x* (t), v* (r) at time r given that it was at x , v at time t may be maximum at 
x* (r) = <x* (r) >,v* (r) = <v* (t) > but it has an ever increasing width asr-t 
increases. If its width were to increase to 277-/k 0 in the time r 2 < r c , the phase 
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of the wave as viewed by the particle would be effectively random and the 
wave-particle interaction would average to zero for r - t > r 2 . In this case 
the width of the autocorrelation would be limited to r 2 . 

Dispersion in statistically averaged orbits in strong turbulence was first 
considered by Dupree. 11 His work is pioneering in the field and shows that 
dispersion leads to a modification of the diffusion coefficient and to resonance 
broadening of the plasma dielectric function. To our knowledge no one has yet 
explicitly considered the effect of modification of the average orbit on either the 
diffusion tensor or the plasma dielectric tensor. 

We have noted that modification of average orbits by the fluctuations and 
dispersion about the average orbits are both effects of order (|8F|) 2 * We expect 
a priori that both affect resonant strong turbulence diffusion with comparable 
magnitude. 

It is our purpose in this paper to calculate simultaneously the effect on 
diffusion of both average orbit modification and orbit dispersion. In a companion 
paper 12 we consider how these two effects non-linearly modify the plasma di- 
electric function and the growth (damping) rate of waves. 

Our approach is via the statistical orbits perturbation technique proposed 
ahd initially developed by Dupree 11 and expanded by Weinstock 13 . As much as 
possible, however, we work explicitly with the dynamic trajectories of the 
particles in phase space rather than resorting exclusively to the operator 
formalism. 
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II. THE DIFFUSION EQUATION FOR STRONG TURBULENCE 


Following Weinstock 13 , let us write Eq. (3) in integral form 


SfU,v, t) - U*C*‘t*)SfC&.)C,6o)- jdvU^SF^z ) . ^1^(7) 


tc 


dv 


Here S f (x , v, t ) is the initial value of S f and the propagator U A ( t , 1 0 ) satisfies 


tr- zUa* (<& **£)■ 


Wa 

T l'*' T ° LJ Zf-W c 5V'/ 
(J4 (£«) £o) - I 

We next make the weak coupling approximation 1 1 • 13 • 14 : 


= 0 


(8) 


Ufl(t,to)2<uct,tc)/>= {&plQ<*cto)-z]-z +[ytio)-y]'jzfy 

(9) 

In Eq. 9 x * and v* are solutions of the characteristic equations 

— * v* PL**M (10) 

<tv d t L ~ y 

again with the boundary conditions x* (t) = x , v * (t) =▼ . The vector U is a 
phase space vector with the six components {x it v i } and has been introduced 
for notational convenience. Note the appearance of F in Eqs. (10) as distinguished 

from <F> in Eqs. 6. Successive higher order corrections to Eq. (9) are in the 
ratio k Q | SF | r c 2 , considered to be a small but finite quantity in the weak coupling 
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approximation.. 13 A straightforward technique for generating higher order 
corrections to Eq. (9) is presented by Birmingham and Bomatici. 14 

The weak coupling approximation is a perturbation scheme for solving 
Eq. (3). It amounts to an iteration on the source term <8F • B S f / Bv>. Elements 
of the convective term 8F- B8 f /Bv are retained in lowest order, for only when 
one includes the influence of fluctuations on particle orbits are the effects of 
strong turbulence manifest. The fundamental strong turbulence solution for 8 f 
represented by Eqs. (7) and (9) thus embodies the fundamental weak turbulence 
solution (i.e., the solution to Eq. (3) with Cl ( 8 2 ) terms discarded) plus select con- 
tributions from all higher order (mode coupling) solutions of the weak turbu- 
lence analysis. In this respect the approximation made here (and also made in 
the work of Dupree 11 and Weinstock 13 ) is similar to those made in other 
theories ( 15 ' 22) of strong plas ma turbulence . 

Having made the approximation 9, we plug the expression for 8 f into Eq. 

(2) and obtain 


r l r ‘ £<& + <E> =■ ~ (<?£' j| fa&tyW&X,** P 


(ii) 


+ 3L . Jdt Lt,v)/S 


to 

We now make the additional standard simplifications: 

(1) We neglect the term in Eq. (11) associated with the initial value 

8 f(x, y , 1 0 ). This term may persist for long times in the phase space 
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description of the plasma, and on this level its omission is unjustified. 
However, when velocity moments of the kinetic equation are taken (and 
we envisage our resultant kinetic equation to be so used), the initial 
value contribution phase mixes to zero in a time interval A t. The 
magnitude of At depends on the nature of 8 f ( x , v, t 0 ). If u, the 
characteristic width of S f (t 0 ) in velocity space is much greater than\ 
the average velocity v 0 associated with S f (t 0 ) , then At = (k 0 u) -1 
where k 0 is the typical wave number of the turbulent spectrum. 

(2) We assume that < f> changes sufficiently slowly with x, v, and t that on 
the right side of Eq. (11) B<f > /3v may be removed from both the <U> 
operation and also from the time integration. Specifically the x and v 
scales of <f> are assumed longer than the incremental change in phase 
space position of a resonant particle during the time r* characteristic 
of the decay of elements of the autocorrelation tensor < S F (x, v , t) 

< U(t, T )> SF (x.vit)> and the time scale of change of <f> is itself 
slower than r*. (We expect that for strong turbulence t* depends on 
the amplitude of the fluctuations in accordance with the discussion in 
the Introduction.) These are standard assumptions 2 3 made in the 
statistical mechanics of stochastic processes. 

With these approximations Eq. (11) again becomes the diffusion equation (4). 

The diffusion tensor is, however, now 
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t-te 




l-t 


uo 


( 12 ) 


" J d%<S£(i,i,t)SP[e*Ct-c), rfi-t), t-t] 7 


The latter form of Eq. (12) follows from the facts that (a) to lowest order in 8F, 
<U (t» t - r) > = u (t, t - r) 11,13 and (b) U (t, t - t) 8F (x, v, t - t) = SF [x* (t -r), 
v* (t - r), t - r] . 14 

If t — 1 0 »r* (Indeed, the diffusion equation is only valid on a time scale 
much longer than r*) the upper limit of integration in Eq. (12) may be extended 
to infinity 

sy 

r ^dr (SF (13) 

0 

Note that the diffusion tensor Eq. (13) differs from the weak turbulence diffusion 
tensor Eq. (5) in the fact that the correlation in (13) is evaluated along the per- 
turbed orbit of the particle. If the turbulence is approximately homogeneous and 
stationary, as we shall assume, the dependence of D in Eqs. (12) and (13) on x 
and t disappears. (The dependence of D on t - t 0 for t - t„ < r* remains 
however). We now attempt to evaluate <u (t, t - t)> in terms of statistical 
properties of the turbulence. 
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in. EXPLICIT EVALUATION OF D 

To abbreviate this calculation, let us restrict ourselves to the case of 
electrostatic turbulence SF = qSE( x , t)/m; the extension to include a turbulent 
electromagnetic field is direct. We further follow Weinstock 13,24 by making a 
cumulant expansion of <U (t, t - r)>, dropping cumulants higher than second 
order: 


c (etp[x\t-z) (e/p A x*(t-zj - 

(14) 

The terms in Eq. (9) involving the velocity increment (v* - v) do not appear in 
(14) because of our assumption of electrostatic turbulence (SF independent of v). 

We plug (14) into (13), expand the SF’ s in Fourier series, and ensemble 
average assuming that the initial phase of each k-mode varies randomly from 
one realization to the next: 


0 ^ 

D - C civ iL k { I ff & I uj l r 


(15) 


exp|-U • ^ i 


< | SFfc | 2> is the ensemble average square amplitude of the Fourier mode with 
wave vector k , k is the unit vector along k, and the imaginary parts of the eigen- 
frequencies are neglected so that the >s are real. The term ik • <Ax* (t - r)> 


10 



contains the effect of the turbulence on the ensemble average orbit. In previous 
work the turbulent corrections to <Ax*> have been neglected so that in an un- 
magnetized plasma <Ax*> was just - vr. The termkk: [<Ax* Ax*> - <Ax*><Ax*>] 
represents the statistical dispersion in the trajectories of the particle 
passing through x, v at time t. 

One can relate Ax* (t - r) implicitly to the fluctuations by formally inte- 
grating Eqs. (10). We here further simplify the problem by assuming that there 
exists no zeroth order magnetic field 

* r v 

ax (t-r) --vv •/- J c(v , T , 3F r £x*(t-ctt'l t-itt'J 

o 

(16) 

A O 

- -vr t J at 1 (t-t) j 

o 

It follows that 


-i)7 = - vr f j dr ' (T- V ) ~rj) - r'^J ' 

C 

r r .t, 
sr -v r r J dr, cfc, 


(17) 


dXfJL-t,) ~ L 


■v 




s - r 1- J r (T- t,; j dx x (t,- rj T 
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in the approximation where terms higher than second order in SF are 


neglected. 

In Appendix A we show that in the random phase approximation T can be 
transformed in such a way that 

, c tr nVl 

= - v* - cic, (t-r,) j 

o 0 (18) 

• SF[x 

For stationary turbulence the correlation in Eq. (18) depends only on | r x - r 2 \ ; 
hence one of the integrations can be performed 


{a A t -1)^ J. ^ jdz (x -z)\sp [x. t ) -SF [x%tX t-zj) 

o ' 

V 

~' V ~ T ~ rjv fa -$Z 


(19) 


The latter approximation follows because the correlation is peaked at z = 0. 
Noting Eq. (12), we finally conclude that 


^ r) j (20) 

where Tr {D} is the trace of the matrix of D. A result similar to Eq. (20) has 
been obtained in one dimension by Weinstock 25 using his operator formalism. 
Our result follows directly from approximating the characteristic trajectory 
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Eq. (16). Note from Eq. (20) that the deviation of the average orbit from a straight 
line may be attributed to an effective fluctuating force (per unit mass) of magnitude 



Next let us calculate the dispersion contribution to Eq. (15). Since 
<Ax* (t - t) > = - vr + C> ( SF 2 ), the dispersion tensor is 


I = <(ax*(V v) a - y^.x \t -rp \*£(t -rj> 


•r 0 r (2i) 

55 1 ^ J ^ ^ \?£ f r ») / t-r-f r, ] c>F + la), t -f rf&jy 


o 0 

Using the stationarity assumption on Eq. (21), we obtain 


I - J c/r, (./ J dtj v z \£F[xh),t'] SF[x*(r- +z£) } v-c, tt£p (22) 

® o 

In Appendix B we show that 


1 ? 


(23) 


Thus the dispersion factor in Eq. (15) is exp - k k: p (v, t) r 3 / 3. Dupree 1 1 
and Weinstock find a similar dependence for the dispersion factor although it 
is the asymptotic D (v,°°) which enters the expressions of these authors. With 
D a function of r the integral in Eq. (15) converges more rapidly than does the 
corresponding integral obtained by Dupree and Weinstock. Dupree's D is also 
by assumption independent of v . 
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We now gather together the pieces, Eqs. (15), (20), and (23) to obtain 



By comparing Eqs. (12) and (13) we see that the non-asymptotic D (v, r) is given 
by the right side of Eq. (24) with r replacing co as the upper limit of integration 


(v,t) J 4t 2 k k ^/df^ | yexp--t[(au £-¥)?'-& jc 77- { § (v, 


(25) 


7 5 

exp - kk ■ 

3 ^ 


Note that Eq. (24) properly reduces to the quasi-linear result 



<*■) - Tl 


Z £ k 


*> 




k v 


(26) 


when the D -dependent exponential terms are suppressed. 

In the strongly turbulent regime, however, we expect that the width of the 
r -integration in Eq. (24) is influenced by these fluctuation-dependent exponential 
terms. The asymptotic D(v, o°) thus requires knowledge of D (v, r), Eq. (25), to 
compute it when the turbulence is strong. Fortunately D (v, r) occurs in Eq. 

(24) in terms which drop off abruptly with r. This fact allows us to use the 
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small r expansion of D (v, r) in Eq. (24). From Eq. (25) we have in the small r 
limit 


£ ‘ 'i) '£ ( tOg - k • v ) J ' t 


(27) 



J) (i/ # r ) = 2 £ \ /% i ? l r ' 


u.U 

4?- 



the term proportional to r 2 vanishing on k -summation. 

We now plug the result Eq. (27) into Eq. (24) and obtain 

fic2 fiFf/exp ^Uuk-k-v) t 


t; 

4 


5'Z ^ '(iSFijj z )( i0 k ' - £ ‘ iOj (28 ) 


.4 /> 


<?>p- — £ * : ^ ^ - 

d t ? / 


Equation (28) is the main result of this theory. 

Let us now look at the time scales on which the integrand of Eq. (28) would 
converge if the three exponential factors occurred individually rather than in 
combination. When contributions from different k’s are summed, the factor 
exp -i(" k - k • v)r produces the usual weak turbulence convergence. This con- 
vergence is due to phase mixing of waves and occurs on the time scale r c = 
(Aw *) -1 • If all waves have approximately the same frequency (as is the case 
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for plasma oscillations, for example), then r c ~ (Ak |y | ) _1 , where Ak is the 
wave number spread in the waves propagating in the direction of v . 

The factor 


Op - 


-~1 

L K k i ?. 
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k P : \£F(xii) Sri*/ 



produces the convergence due to the statistical spread in particle orbits. This 
convergence occurs on the time scale r 2 = (m/k 0 q SE rms ) !4 . r 2 is the time 
it takes the dispersion in the orbits of a particle to extend over one typical wave- 
length 2n/k 0 of the fluctuations given that the orbits start at the same phase 
space point in all realizations. (In Dupree’s notation r 2 is r tr .) 

Finally, the term 



produces the convergence associated with the acceleration of the ensemble aver- 
age orbit on the time scale 



t'n 






In estimating t x we have taken the typical wave frequency seen by a resonant 
particle to be of the order Aw*, the frequency spread observed by the particle. 
The complicated parametric dependence of r x reflects the fact that the effective 
fluctuating force 3/Sv • JD behaves as 
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~ 3 a ,* 

C AK ^ id 


SL Ml 

hi 


, r^ r 6' ■ 1 


Since r 2 depends on (SE rms ) -i4 , we expect that for very large fluctuations 
r 2 « r c . Unfortunately the weak coupling approximation 11, 13 breaks down for 
r 2 < r c and hence Eq. (28) is invalid in this regime. 26 

Equation (28) can, however, be examined for r c « r 2 and the effects of 
average orbit distortion and orbit dispersion assessed. As previously men- 
tioned, these orbit effects act to modify the coherence time between the fluctua- 
tions and resonant particles. Note that when r c « r 2 and Ak ~ k 0 , r t is a 
shorter time interval than r 2 . Thus average orbit distortion is a more important 
effect than orbit dispersion in determining corrections to D in this regime. 

Ak ~ k Q implies that the spectrum is broad-band, as would result, for example, 
from a beam, diffuse in energy, penetrating and giving rise to instability in a 
background plasma. 

Figure 1 illustrates the effect of average orbit distortion on the diffusion 
coefficient for a one dimensional plasma. The turbulence is here characterized 
by a spectrum 

S I c n | a. ^ j [~ (A- A>) 1 ' a [ At A> ) 

<lw *. I / 04 f 


wherein k 0 /Ak has been chosen to be 4. We have also assumed that the frequency 
of the component waves is independent of wavenumber: co k - cl> 0 for k > 0 and 
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oj k = -co 0 for k < 0. Plotted against particle velocity (in units of a> 0 /k 0 ) is the 
non-dimensionalized diffusion coefficient. 

Curve I is the usual weak turbulence coefficient Eq. (26). Curve II includes 
the effect of average orbit distortion Eq. (28) ; dispersion has been neglected. 
(Normalization of the two curves is the same, so that they may be compared in 
amplitude as well as shape.) For the spectrum we have assumed, 

2 -k 1 ISP k il a ‘L UJ V-L'v') l/(/-/'03#) (29) 

A 1 

where 

vv ~ 2 s’Fcx-j 

is the mean square value of the fluctuating force (per unit mass). In obtaining 

i *“<f 2 

Curve II we have had to specify a value for the parameter e = -r- ; 

"o 2/k o 2 - 

we have taken e to be 10 -3 . e measures (in units of the resonant energy 
-o 2 

m — - ) the energy that a particle would gain in one period « 0 1 of the fluctua- 
k o 

tions if it were to freely accelerate in the rms electric field. Our choice of e 
insures that for resonant particles t c < t 2 . 

Note that Curve II has an extremely rapid initial fall-off from its maximum 
value at v = 0.97. (The maxima of the two curves coincide for our choice of 
parameters.) This behavior is attributable to the factor (1-1.03 v) in Eq. (29). 

If this factor were equal to unity — or if dispersion were added — the maximum 
of Curve n would be both lower and more rounded. (We feel that dispersion is 
important only in the immediate vicinity of the maximum.) 
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For Ak = k 0 , the bulk of the resonant particles lie between v = .8 and 

v = 1.3. For particles in this region the effect of average orbit dispersion is 
generally to reduce the magnitude of the diffusion from its weak-turbulence 
value. On the flanks of the resonance, v < .8 and v > 1.3, the diffusion coeffi- 
cient of Curve II is larger than that of Curve I: the turbulence is here evidently 
acting to make non-resonant particles more resonant. 

SUMMARY 

We have derived the strong turbulence diffusion tensor D. D includes the 
effects of fluctuations on both the ensemble averaged plasma orbits and the 
statistical dispersion in these orbits. By using the plasma trajectory equations 
we have calculated the contributions of these two effects explicitly. 

There are two basic time scales in this problem: r c is the autocorrelation 
time of the fluctuations as seen by a resonant particle when fluctuation effects 
on its orbit are ignored and r 2 is the time it takes a resonant particle to diffuse 
(under the influence of the turbulence) one characteristic wavelength (2 7r/k 0 ) of 
the fluctuations. In the weak coupling approximation to strong plasma turbulence 
we consider r c « t 2 . 

Effects of average orbit distortion occur on a r x time scale, where 

/ k 0 \ 2/5 

r i ” VAk/ ( t 2 4t c) 1/S ' Thus when k 0 ~ Ak, r x < r 2 in the weak coupling 
limit and the effects of average orbit distortion dominate those of orbit dispersion. 

In a sample situation studied in this paper, average orbit distortion reduces 
the magnitude of the diffusion for the bulk of resonant particles. The diffusion 
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coefficient for particles on the flanks of the resonance is, however, enhanced 
as the fluctuations here apparently act to make particles more resonant with 
the turbulence. 
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APPENDIX A 


The Correlation T 


Expanding the fluctuating force 8F in a Fourier Series, we express T as 


A A. 


T - (iZ k‘lsF l 'llSFk'\ 

e/f L [k-x % -Zl) r to* [t -tx) f 4>h 

t-k 1 ' X fOJp i-b-li) f 4' J 


(A.l) 


As in Eq. (15) the £ 's are unit vectors in the direction of electrostatic wave 
propagation, | SF^ | is the rms amplitude of the mode with wave vector k , and 
is the phase of the same wave at x = 0 at time t = 0. The 0 k ’s are defined 
so that the k’s all have positive projection on an arbitrarily chosen axis. We 
interpret the ensemble average as an average over the phase angles <p k . If we 
assume that for each k, c£ k varies randomly from one realization to the next, it 
follows that only the terms k 1 = -k contribute in Eq. (A.l). 


T = -i 2. * l$r t . j z fyp i ■[**&- ' 'x.(±-£i)j - u>£ (r 2 -n) v < A * 2 > 


Realize now that, by Eq. (16), 


X*tt -tx) - ) =■ fc -?*■ J ~ &X *C* ~ti ) 


(A.3) 



so that 


Jl k e*p Ck , [ ** it -Vz) - 


(A.4) 


7 , h ' l i- OQ?) 

iSarU) ~ / 


It thus follows that 


i rj . . _ r 7 

? t k F !^<*l-^p-kki t,)-X%-T,) J-WttH-ty (A.5) 


But in the random phase approximation 


{SP [x*(t -r 3 ) ; t-lji j - J'f [&*&*,), i- rJ/ > - 


and 


T , 


(A .6) 


z te / i( & ' [xV^j - -ti)]- Uj k j 


/ cP 


T a -r, 


~ ^ sr[xJi±-^),t-Za] ■ Jr [x*ci r) } -t-ry <a.7> 


<Ax* (t - r)) thus has the form Eq. (18). 


25 



APPENDIX B 


Evaluation of I 


To evaluate Eq. 22, let us first introduce the variables 


S~ Z, 



(B.l) 


By the stationarity hypothesis <8F[x*(r),r) SFtx^T-Tj + t 2 ), r-r x +r 2 ]> 
depends only on | s |. Since r 1 r 2 = (q 2 - s 2 )/4, we need hence only integrate 
Eq. (22) over the region s > 0 (r t > t 2 ). Thus 



4t -S 


<*1 ^r[x*R)/c) Sf 


(B.2) 



But the correlation <SF[x* (t),t] SF[x* (r - s ),r - s] > is strongly peaked at 
s = 0, so that 


I~r c *J 42 \w[xfy,t]sF[x*(t-s) ? T-€s> 

o ~ x 

= |r 3 J ds \SP[x/t) S£ 


(B.3) 


the latter form of (B.3) again following from the stationarity assumption. 

Compare now the integral in (B.3) with the definition of the diffusion tensor 
Eq. (12). They are the same and Eq. (23) follows. 
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